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UNIFORM SPECTRAL AMPLITUDE WINDOWING FOR 
HYPERBOLIC-FREQUENCY MODULATED WAVEFORMS 


INTRODUCTION 

The development of classification algorithms within the Naval Research Laboratory Shallow 
Water Active Classification Project (ONR Project RJ14B87) requires that the response of the 
target be measured over a broad range of frequencies. This must be done with sufficient fidelity 
to resolve the subtle acoustic features of the target. It will allow algorithms to be developed that 
determine not only the class of target present (if one is present), but also target aspect. 

A waveform used in an active sonar system that collects data for a classification algorithm 
must meet several criteria. First, the waveform must be broadband and fill all available bandspace 
offered by the transducer. The waveform shovtld produce a high-resolution range profile of the 
target when matched filtered, so it must be a pulse compression waveform. Moreover, it should 
be Doppler-tolerant, that is, the zero-Doppler matched-filter response to a return should yield 
the range profile of the target even if the taiget is moving. Finally, the waveform should have a 
flat spectrum; this will allow a quick in situ look at the target spectrum without extensive post 
processing to compensate for the spectral coloring of the waveform. It wiU also simplify the design 
of a classifier. 

With these criteria in mind, we consider the hyperbolic frequency-modulated (HFM) waveform 
as a good candidate. In its traditional form, it immediately meets most of these criteria; it is a 
broadband, Doppler-tolerant, pulse-compression waveform [1-5). The waveform does not have a 
flat spectrum, especially below 1 kHz. However, with the correct time window, the spectrum can 
be made esseitially flat. The design of such a window is the subject of this report. 

THE HYPERBOLIC FREQUENCY-MODULATED WAVEFORM 

The time-frequency function of the hyperbolic frequency-modulated waveform is given by 


where /(O) = /i, /(T) = / 2 , P = /a — /i is the bandwidth (in Hz), and T is the time length 
(in seconds) of the waveform. Note that Eq. (1) is a hyperbolic function, hence the name of the 
waveform.* If /i < /a, then the waveform is an upsweep, otherwise, it is referred to as a downsweep. 
Waveforms using this form of time-frequency function possess the same pulse-compression property 

MunNcript if^iDved Fsbniaiy 1,1904. 

*Note that the wavefonn period given by l=f(t) is a linear function, thus, the waveform is also known as a linear 
period modulated (LPM) waveform. This name is commonly iised among those working in high-frequency acoustics. 
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offered by the matched-filter processing of a linear frequency-modulated (LFM) waveform, but they 
are also Doppler-tolerant. That is, the matched-filter response to an echo from a moving point target 
will have a large peak value, even if the return is processed under the assumption that the target 
is stationary. 


The functional form of the waveform with a rectangular window is given by 


x{t) = 


( 


0 


ift€(0,ri, 

otherwise, 


( 2 ) 


where >4 is an arbitrary constant. Prom Eq. (1) the instantaneous phase is 

<P{t) = 27r S{r)dT = ln(/ 2 r - Bt ). (3) 

Figure 1 shows the spectrum of an HFM waveform for /i = 500 Hz, = 700 Hz, and T = 2 
seconds. Note the distinctive slant to the spectrum. This feature is most prominent in practical 
HFM waveforms whose spectral centroid is below 1 kHz. In the next section we derive a time 
window that will remove this feature of the spectrum. 



Fig. 1 - Spectrum of an HFM waveform with /i = 500 Hz, /2 = 700 Hz, 
T = 2 seconds, and a rectangular window 


WINDOW DESIGN FOR A FLAT SPECTRUM 


Consider the case of an HFM waveform whose time-frequency function and phase are given by 
Eqs. (1) and (3), respectively, but whose time function is now given by 


x{t) = 


1 


a(t)eJ^‘> iffe{0,ri, 
0 otherwise, 


(4) 


where o(f) is an arbitrary non-negative, real window. We now show that it is possible to chose a{t) 
such that the waveform’s spectrum will be flat. Here, the term ‘flat’ is used in a qualitative sense, 
sinoe the example given below will show that the spectrum can still have a large amount of ‘ripple.’ 
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Pig. 2 - Time-frequency fimction of a frequency-modulated waveform, and 
the telation«hip between its alope, the dwell time At, and the frequency 
band Af 


Consider the arbitrary time-frequency function shown in Fig 2, denoted by f(t). Also, let 


nit) = dfit)/dt. 


(5) 


Consider a band of frequencies of size A/ centered about f{t), then the approximate ‘dwell time’ 
of the waveform within this band is approximately given by 


At = 


A/ 
ln(t)l ■ 


( 6 ) 


It follows that the approximate spectral energy in the band A/(t) is 


E = a\t)At. (7) 

Since we want the spectrum to be flat, the'amount of spectral energy in all bands of fixed size A/ 
should be constant, irrespective of our choice of t. Thus, from Eqs. (6) and (7) it follows that 


a\t) 


A/ 

mi 


= constant. 


( 8 ) 


Equation (8) implies that our choice of the window requires 


o(0 « y/m\- (9) 

Note that this expression is not specific to an HFM waveform; it is valid for any frequency- 
modulated waveform. However, for the specific case of an HFM waveform, combining Eqs. (1), 
(5), and (9) yields 

^ /2 - Bt/T ’ 

where K is an arbitrary constant whose value depends on the application of the waveform. If the 
waveform is to be used in a peak-power-limited system, it must obey the property 0 < a(t) < 1, 
implying that 
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Fig. 3 - Spectrum of an HFM waveform with /i = 500 Hz, /* = 700 Hz, 
T = 2 aecondt, and a window designed to remove spectral slant 


^ h- Bt/T ^ ^ 

If, however, the waveform is to be used for performing matched-filtering, then the energy normalized 
window can be found by integrating the square of the window given in Eq. (9), setting the result 
equal to 1, and solving for K. The result yields 


o(t) = 


h-BtlT' 


( 12 ) 


Figure 3 shows the spectrum of an HFM waveform with the same frequencies and time length 
as the waveform whose spectrum is given in Fig. 1; however, we have applied the energy normalized 
window given in Eq. (12) and shown in Fig. 4. The spectrum is quite similar to that of a linear 
frequency-modulated (LFM) waveform. 


REDUCTION OF SPECTRAL RIPPLE; COMPOSITE WINDOWING 

Although the ‘slant’ to the spectrum of an HFM waveform can be removed by applying the 
window derived in the previous section, Fig. 3 shows that a large eunount of spectral ripple cam still 
be present. The ripple is also known as Gibb’s phenomenon amd is due to the presence of the jump 
discontinuities in the window at t = 0 and t = T [6,7]. 


To reduce the spectral ripple, we propose the use of an additionail window to maJce the wave¬ 
form’s functional definition continuous. In this report we use the tapered-cosine window given 




t € [0,ar), 
t€lar,(l-a)T), 
t€l(l-a)T,r], 
otherwise, 


(13) 


k.... 
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Fig. 4 - Envelope of an HFM waveform with /i = 500 Hz, /j = 700 Hz, 

T = 2 aeconda, and a window designed to remove spectral slant 

where a is a parameter that determines the effective width of the window, and 0 < a < 0.5. Use 
of the window in Eq. (13) implies that the waveform’s functional form is now given by 

x(t) - ! if f e [0,T], , . 

~ \ 0 otherwise, ^ ^ 

where a(t) is given by Bq. (10). There are, of course, other possible choices for w(t), but the form in 
Eq. (13) was chosen because it is a common window type. It is also once continuously differentiable 
with a bounded second derivative, which implies that the waveform’s spectrum obeys the property 

X{w) = ~ ^ 00, (15) 

where is the spectral centroid. This property implies that the ‘spectral skirt’ will be lower than 
that of the spectra for the waveforms given by Eqs. (2) and (4). 

Figure 5 shows the spectnun of an HFM waveform with the same frequencies and time length 
as in Fig. 1. However, it has the functional form given in Eq. (14), for which the composite window 
a{t)w{t) is shown in Fig. 6, with a = 0.1. The waveform was energy normalized numerically. 
Clearly, the magnitude of the spectral ripple has been reduced. 

If the waveform is to be used in a peak-power-limited system, then we must find the constant 
K such that 0 < w{t)a(t) < 1, where o(t) is given by Eq. (10). To do this, we first determine where 
the peak value of the window should occur. Some thought reveals that for an upsweep the peak 
should occur in the interval t € [(1 — a)T,T\, and for a downsweep the peak should occur in the 
interval t € [0, oT]. Furthermore, the derivative of the composite window must be equal to zero at 
the peak. Hence, differentiating a(t)ti;(t), and setting the result equal to zero leads to the following 
equation whose solution gives the location of the peak if the waveform is an upsweep (/i < /a): 

Bo|l-cos|^^^l-|;^JJY)sinj^^^l-|;jJ =0 for t € ((1-a)r,ri. (16) 
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frequency (Hz) 

Fig. 5 - Spectrum of an HFM waveform with /i = 500 Hz, ft = 700 Hz, 
T = 2 seconds, and a composite window with a = 0.1 



Fig. 6 - Envelope of an HFM waveform with fi — 500 Hz, ft = 700 Hz, 
T = 2 seconds, and a oompomte window with a = 0.1 
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By a similar procedure, it is possible to derive the following equation whose solution gives the 
location of the peak if the waveform is a downsweep {/i > / 2 ): 

Ba [^1 - cos +^(j 2 - y) s® (^) =0 for t€ [0, aT] (17) 

Denoting the solution to either Eq. (16) or (17) as tp, it follows from (10) that K = f 2 — BtpJT. 

Obviously the equations in (16) and (17) are nonlinear and must be solved numericeilly. How¬ 
ever, several numerical algorithms can be used to find the solution, such as the bisection method, 
secant method, regula falsa, and Newton-Raphson. Luckily, some relief from this problem can be 
found by using mathematical software packages such as Mathematica [8] or Maple (9) that have 
these numericsd algorithms preprogrammed as functions. 

BANDWIDTH COMPENSATION 

Use of the tapered-cosine window will trim the leading and trailing edges of the envelope, and 
consequently remove those spectral components associated with the beginning and end of the pulse. 
The result is a reduction of the true bandwidth of the waveform; this can be seen by comparing the 
spectra shown in Figs. 3 and 5. In this case, the window effectively removed the leading and trailing 
5% of the energy in the window, hence, the spectrum in Fig. 5 is approximately 10% narrower than 
the spectnun shown in Fig. 3. 

To design a waveform with the desired bandwidth B, we must adjust the values of the fre¬ 
quencies /i and /z in Eq. (1). In this case, we define the ‘design frequencies’ f\ and /z such that 
B — f 2 — fi, and now refer to /i and /z as the ‘parameter frequencies.’ The design frequencies 
are known. The parameter frequencies are used in Eqs. (1) and (3) and must be derived from the 
design frequencies. 


To derive design equations for the parameter frequencies we make the following observation: 
the tapered-cosine window as given in Eq. (13) effectively trims the leading and trailing 100(a/2)% 
of the energy from the rectangular window, so the bandwidth is approximately 100a% of that 
for the rectangularly windowed HFM. To compensate for this effect, we set fi = f{aT/2) and 
/a = /((I — Q!/2)T). Thus, from Eq. (1) and the definition of the bandwidth B, it is possible to 
show that 



The simultaneous equations in (18) can be easily solved for the unknown values /i and /z by using 
Gaussian elimination, or the method of determinants. The solution is 


2(1-Q) 

{2-a)/h-a/f2' 

2(1-g) 

{2-a)/f2-a/h 


(19) 
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Pig. 7 - Spectrum of an HFM waveform with fi = 500 Hz, /j = 700 Hz, 
T = 2 seconds, and a composite window with a = 0.1 


We continue with the examples used in the previous sections based on an HFM waveform with 
/i = 500 Hz, /2 = 700 Hz, and a = 0.1. Substituting these values of the design frequencies into the 
equations in (19) yields fi = 492.19 Hz and = 715.91 Hz. Figure 7 shows the spectrum of a unit 
energy waveform using these parameter frequencies. It can be seen that it has the same bandwidth 
as the waveform that uses a rectangular window, whose spectrum is shown Fig. 3. 


The autocorrelation function of the waveform x(f) is defined as 


fl(r) = \f_ 


‘OO 

00 


x(t)x*(t — T)dt 


( 20 ) 


Fig. 8 shows the autocorrelation fimctions of two of the HFM waveforms given in the examples. The 
dashed line is the autocorrelation of the rectangularly windowed HFM whose spectrum is shown 
in Fig. 1. The solid line is the autocorrelation function of the HFM presented in this section, to 
which we have applied composite windowing and bandwidth compensation. Note that the shape 
of the main peak and the location of the sidelobes 6ire essentially the same for both waveforms. 
This demonstrates that the use of the composite window and beuidwidth compensation can have 
negligible effects on the resolution of the waveform. 


CONCLUSIONS 


We have presented a method of designing HFM waveforms that have essentially flat spectra 
possessing a small amount of ripple. This was accomplished by selecting an appropriate time 
window. The form of the window o(f) given in Eq. (10) was also suggested by Kroszczynski in 
Ref. 4, but its derivation is contained in a reference cited in Ref. 10. It is not known if the derivation 
of the window is similar to the approach presented here, since the reference is not available. 
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delay (s) 


Fig. 8 - Autocorrelation functions of two HPM waveforms. Dashed line is 
the autocorrelation of an HFM waveform with ft = 500 Hz, /j = 700 Hz, 
T = 2 seconds, and a rectsmgular window. Solid li^e is an HFM with the same 
design frequencies and time length but with a composite window (a = 0.1) and 
bandwidth compensation. 
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